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ABSTRACT 

We explore the chemical structure of a disk that contains a large central gap of R ~ 45 AU, as is 
commonly seen in transitional disk systems. In our chemical model of a disk with a cleared inner 
void, the midplane becomes revealed to the central star so that it is directly irradiated. The midplane 
material at the truncation radius is permissive to reprocessed optical heating radiation, but opaque to 
the photo-dissociating ultraviolet, creating an environment abundant in gas-phase molecules. Thus the 
disk midplane, which would otherwise for a full disk be dominated by near complete heavy element 
freeze-out, should become observable in molecular emission. If this prediction is correct this has 
exciting prospects for observations with the Atacama Large Millimeter/Submillimeter Array (ALMA), 
as the inner transition region should thus be readily detected and resolved, especially using high-J 
rotational transitions excited in the high density midplane gas. Therefore such observations will 
potentially provide us with a direct probe of the physics and chemistry at this actively evolving 
interface. 

Subject headings: accretion, accretion disks — astrochemistry — circumstellar matter — stars: pre- 
main sequence 



1. INTRODUCTION 

Over the last thirty years, our knowledge of pre- 
main sequence evolution has undergone significant ad- 
vances. Observations of the full spectral energy dis- 
tribution (SED) from disks around young stars have 
shown that not all are alike, with the striking discov- 
ery of a subset of disks with optically thin inner "holes" 
devoid of small grains, s urrounded by optically thick 



Strom et al. 1989; Calvet et al. 2005; 



Espaillat et al.ll2010f ). This inner void 



outer disks (e.g 
iCieza et al.l [201(1 
was later confirm ed by r esolved sub-mm i nterferometry 
(e.g. iPietu et al.l[2 006: B rown et al.l 120081 : iHughes et al.l 
120091 : [Andrews et al.l l200l l2011f ) and was interpreted 
as an intermediate stage betwee n primordial and debris 
disks, coined "transition disks" ([Strom et al.lll989l ). 

In this work we adopt the physical definition of "tran- 
sition disk" as an optically thick disk truncated within 
some inner radius within which there has either been 
substantial grain growth or removal. Such objects 
are of particular interest, as the presence of a gap 
has been attributed to c l earing by young protoplanets 
fe.g.lSkrutskie et al . [[199(1 [Brvden et al.[[1999[ : rRice et al.[ 
[2003f) or tidal inter actions with young stellar companions 
( e.g.^utrev et al.[[l99l). 

[Dullemond et al.[ ()2001[ ) investigated the physical 
structure of Herbig Ae/Be disks that possess large inner 
radii similar to those seen in low mass transition disks. 
In this work it was found that at the inner rim, the nor- 
mal incidence angle of the stellar irradiation causes the 
disk rim to be much hotter than it would otherwise be 
for a classical flared disk, where the radiation arrives at 
a glancing angle. Further observational evidence for the 
directly irradia ted wall was later see n both in studies of 
the SED (e.g. [Espaillat et al.[[2007[ ) as weh as through 



scattered light imaging ([Brown et al.[[2008[ ). 

ALMA will readily resolve such inner gaps and voids 
in the dust disk (e.g. iWolf et al.l[200llWolf D^Angeld 
[2005[ ): however, here we seek to make predictions regard- 
ing the significant gas reservoir in transition disks. In 
particular, one interesting aspect of disk chemistry that 
has yet to be explored is the potential that this warm UV 
irradiated inner rim should have unique chemical proper- 
ties. This has important implications since it is possible 
that the physical and chemical conditions of the previ- 
ously hidden, but now exposed, midplane will be revealed 
and potentially detectable with high spatial resolution 
observations, such as those anticipated by ALMA. Here 
we present a chemical model of a protoplanetary disk 
with a large inner gap si milar to those s een in classical 
trans ition disks (e.g. Hughes et al.[[2QQ9[ : [Espaillat et al.[ 
[2010[ ) as a prospective study of observability with ALMA. 

2. CHEMICAL MODELING 

2.1. Disk Framework 

For our model we ado p t the comprehensive disk struc- 
ture from ^ Sauter et al.[ ([200 9*). originally purposed for 
CB26 and constrained by a large set of observations: the 
SED from nanometer to millimeter wavelengths, near- 
infrared scattered light images, and resolv ed millimeter 
image s. For the central star, the model of [Sauter et al.[ 
(|2009[ ) assumes standard T Tauri values of M = 0.5 M© 
and L = 0.9 L©. The disk density profile is given by: 
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Usin^ this model a nd temperatu r es der ived with MC3D 
([Wolf et alJll999[ ). iSauter et all (j2QQ9[ ) successfully re- 
produced the full set of observations with the following 
best-fit model parameters: ro = 100 AU, ho = 10 AU, 

a = 2.2, P = 1.4, i^outer = 200 AU and i^inner = 45 

AU. The measured size of the hole is typical of that 
seen in resolved millimeter observations and SED mod- 
eling of transition disks, which range from a few AU to 
upwa rds of 70 AU (jAndrews et al.ll201ll : lEspaillat et al.l 
120071 ). The inner disk irradiation substantially increases 
the disk temperature at the wall, reaching ~ 50 K, as 
compared to a typical midplane temperature of ~ 15 
K as seen at larger radii. We note that this model is 
taken as a "snapshot" and any subsequent disk physical 
evolution is beyond the scope of this work, though the 
chemistry of a variety of disk geometries will be explored 
in an upcoming paper. 

For the gas density we assume a standard ISM gas- 
to-dust mass ratio of /g = 100 and that the dust and 
gas are co-spatial. The gas temperature is taken to be 
equal to Tdnst, which hold s true near the dense mid- 
plane (*Jonk heid et al.|[2'004f ). We furthermore note Tdust 
was originally derived assuming passive heating by the 
central star, with a ccretion heating treated as negligible 
(|Sauter et al.ll200"9l ). The effects of inclusion of accretion 
heating are discussed in Section |4]TJ The model physical 
structure is shown in Fig. [H 

For the opacities we adopt IWeingartner fc Draind 
(|2001f ) for a blend of astronomical silicates and carbona- 
ceous grains with Ry = 5.5. We also assume an unset- 
tled disk and uniform dust composition for simplicity. 
Deviations from these assumptions, such as the affect 
of vertical settling o f sm all grains and grain-growth, are 
discussed in Section l4?2l 

2.2. Radiation Field 

Both the ultraviolet and X-ray radiation field are be- 
lieved to be dominant factors driv i ng the chemistry in 
disks (jGlassgold et al.l I1997L 120041 : Ivan Zadelhoff et al.l 
[2003h . While low mass young stars peak in the opti- 
cal regime with negligible chromospheric contribution to 
the FUV, the accretion shock at the central star provides 
a significant source o f both FUV and soft X-ray flux (e.g. 
iGiinther et al.l [20071 ). We neglect external contribution 
from the ISRF since the stellar FUV will dominate ex- 
ternal sources by o rders of magnitude in the inner disk 
(jBergin et al.ll2003[ ). though its import may lie in under- 
standing the outer disk chemistry (Obe rg et al.l |201Q| ). 
Due to the geometry of the star-disk system it is then 
neces sary to properly treat the radiative transfer into the 
disk (IWillacv fc Langeij 120001: Ivan Zadelhoff et al.l[200l 
iFogel et al.l 1201 If ). This is especially true in transition 
disks, where voids and gaps can allow photons to prop- 
agate more freely into the outer disk material. Possible 
implications of the presence of a small amount of unde- 
tected dust within the gap are discussed in Section l4?3l 

2.2.1. Continuum Radiative Transfer 

We assume the FUV continuum opacity is dominated 
by the dust, and consequently the UV fleld is depen- 
dent on settling of small grains, opacities assumed, and 
disk geomet ry. Using the measured FUV spectrum of 
TW Hydra (|Herczeg et al...2002. .2004 ) binned down to 
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Fig. 1. — Left: Plot of the transition disk model with the cen- 
tral star at the origin. Panels: a) gas density b) temperature c) 
FUV radiation field at 1600 A. Right: Vertical cuts of the same 
quantities taken at R = 45, 50 and 100 AU. 

nine discrete wavelengths between 950-2000 A, we cal- 
culate the continu um radiative trans f er int o the disk us- 
ing the method of iBethell fc BerginI ()2011l ). implement- 
ing the opacities described in Section 12.11 The choice 
of nine wavelengths was motivated by the need to suffi- 
ciently capture the shape of the FUV continuum, to be 
insensitive to individual weak emission lines, and to also 
remain computationally efficient. 

2.2.2. Lyman-a Radiative Transfer 

While many individual lines in the observed TW Hy- 
dra spectrum are weak, the Lym an- a line alone car- 
ries ry 85% of the total FUV flux (jHerczeg et al.l 120041 : 
iBergin et al.ll2003"[ ). This line is furthermore expected to 
play a significant chemical role as a number of molec- 
ular species have photodissociation cross sections near 
1216A (|Fogel et al.l 1201 If ). In addition to dust scatter- 
ing, Lyman -g undergoes isotropic scattering off hydro- 
gen atoms (jBethell fc BerginI |201 If ), which in principle 
requires an iterative and computationally expensive cal- 
culation between the radiation field and the chemistry. 

In this work, we incorporate the effects of Lyman- 
a using an approxim ate treatment motivated by 
IBethell 



BerginI ()2011l ). The procedure is illustrated 



in Fig. f2l If one assumes a priori that there is an opti- 
cally thick layer of atomic hydrogen on the disk surface 
located at rpuv = 1 as defined by the dust (Fig. |2l Re- 
gion II), the Lyman-a photons will first encounter this 
isotropic H-scattering layer. A fraction of the radiation 
will be lost to space (< 50%, Region I: wavy arrows), 
with the surviving Lyman-a photons diffusively propa- 
gating through the H-layer. Below the atomic layer, the 
hydrogen is predominantly H2 and the Lyman-a pho- 
tons proceed as continuum photons scattering off only 
the dust grains (Fig. f2l Region III). Thus, at the base of 
the hydrogen scattering layer (II^III), the photons effec- 
tively form a layer of isotropically emitting point sources, 
which together create a planar source of radiation that 
shines at a nearly normal angle to the midplane. There- 
fore even though a significant fraction of Lyman-a pho- 
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Fig. 2. — Lya schematic illustration. (I.) "Free-streaming" re- 
gion above Tradial ~ 1 surfacc; half of the Lya radiation escapes the 
disk. (II.) H-scattering layer with a vertical depth of Tvertical ~ 1- 
Remaining photons diffuse downward until they escape the atomic 
region. (III.) Region where H is predominantly H2. Dust now be- 
comes the dominant source of opacity to Lya photons, and the Lya 
radiation now behaves as vertically attenuated continuum photons, 
which significantly enhances its penetrating power. The (I.) — )► 
(II.) interface illustrates the r ~ 1 surface for continuum photons, 
whereas the (II.) — )► (III.) interface is effectively the r ~ 1 surface 
for Lyck; photons. 

tons are lost to space, the remaining photons wih have 
greater vertical penetration power into the disk and can 
propagate many AU deeper than the UV continuum ra- 
diation. On the front edge of the disk, at the dense inner 
rim, Lyman-a is similarly "stopped" as was seen for the 
continuum photons, (cf. Fig. [T] (c)). We discuss the 
implications of an inner void that is not empty of gas on 
Lyman-a transfer in Section [4l4l 

2.3. Reaction Network 

Combining the model detailed in Section 12.11 and the 
radiation field in S ection [2.2[ t he re sultant chemistry 
is calculated using iFogel et al.l ()2Qllf )'s comprehensive 
disk chemical model, based on the Ohio State Univer- 
sity Astrophysical Chemistry Group's gas-phase network 
(jSmith et al ll2QQ4 ). The reaction types include photo- 
desorption, photo-dissociation, freeze-out, grain surface 
reactions, ion and electron reactions, cosmic-ray and stel- 
lar X-ray ionization, and radiative reactions. In total, 
the network encompasses 5910 reactions and 639 react- 
ing species, including some time-dependent reactions, 
encompassing the main species of astrochemical im por- 
tance, and described in detail in lFogel et al.l ()2Qllf ). 

The model initially assumes uniform niolecu lar cloud 
chemical abundances lAikawa fc HerbstI () 19991 ) and fol- 
lows the chemical evolution for 3 Myr. In this work we 
take the abundances at 1 Myr, which is long enough such 
that the chemistry has "relaxed" but not so long that the 
disk would likely have physically evolved away from this 
state. 

Furthermore, we assume a typical integrated T Tauri 
star X-ray luminosity of 10^^ erg s~^ and a thermal X- 
ray spectrum between 1-10 keV (jClass gold et al.l [ 19971 
and references therein). Th e model of[F ogel et al. (2011) 
incorporates the method of lAikawa fc Herbst (2001 ) for 
X-ray propagation. For cosmic rays we adopt a typ- 
ical cosmic ray ionization rate of 1.3 x 10~^^ s~^ 
per H, with an attenuation column of 96 g cm~^ 



(|Umebavashi fc NakariQlll981f ). 

3. RESULTS 

3.1. Disk Chemistry 

The resultant chemical abundances (relative to nu = 
nni + 2nH2 ) observed gas-phase species are shown 

in Fig. [31 Molecules such as CO, H2CO, and N2H+ show 
an enhanced gas-phase abundance at the wall. For com- 
parison, in an untruncated disk model, all neutral species 
plotted in Fig. [3] would be otherwise frozen onto grains 
at the midplane at Rwaii = 45 AU. The two ions shown, 
N2H+ and HCO^, would also not be present in an un- 
truncated disk, as their chemical precursors, N2 and CO, 
would be frozen out, inhibiting the formation of these two 
species. Consequently the gas-phase enhancement at the 
truncation radius shown for all species plotted here, be- 
sides H2O, would not exist if it were not for the large 
inner gap. 

Water, with a freeze-out temperature of ~ 100 K, is 
frozen onto grains at the ~ 50 K transition region. One 
could however envision a disk with a wall closer to the 
central star and thus warmer, such that water would sub- 
limate from grains and be observable. Thus the specific 
species present are not necessarily the most important 
result, but that the star can efficiently heat the wall, 
warming it above the sublimation temperature of a vari- 
ety of species, potentially allowing us to observationally 
probe disk physics as well as evolutionary state. There- 
fore the existence of a cleared inner gap should produce 
unique chemical features in the outer disk not present in 
full classical disks. 

3.2. Observables 

If such an enhancement is present it is of interest to 
determine observability. While previous observ ations 
reveal a diverse chemistry (e.g. Dutre y et all 119971 : 
lOberg et al.l l201Cil ). the sensitivity and resolving power 
of ALMA is required to fully understand the detailed 
structure of these systems. 

We adopt a disk inclination of 60° and calculate the 
resulting emission for rotational transitions of ^^CO and 
0-H2CO as would be seen at a distance of 140 pc us- 
ing the non-LTE line radia tion transfer code, LIME 
(|Brinch fc Hogerheiidel[20Toh . These species have been 
chosen as both are commonly observed towards disks 
aroun d low mass T Tauri stars (e.g. lOberg et al.ll2010l . 
I2OIII ). The densities reached at the frontally illuminated 
midplane preferentially excite high-J rotational transi- 
tions, and therefore transition disks should uniquely ex- 
hibit high-J bright molecular rings at the wall when re- 
solved by sub-mm observations with ALMA. Using this 
calculated line emission, we then use the SIMDATA pack- 
age in CASA to calculate the ALMA visibilities and re- 
construct the image from the UV coverage of the full 
ALMA array. Fig. |4] shows the simulated observations 
for these lines as seen by ALMA for an antenna con- 
figuration with ~ 0.1^' resolution at 672 GHz, chosen 
to provide adequate resolution and sensitivity assuming 
typical thermal noise. 

4. FURTHER CONSIDERATIONS 

In this work we propose the presence of gas-phase 
molecules at the inner edge of transition disks. This ef- 
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Fig. 3. — Chemical model results, plotted as abundance relative to the total number of hydrogen atoms. Shown are common species of 
astrophysical interest: (a) CO, (b) HCO+, (c) H2CO, (d) H2O, (e) N2H+, (f) HCN. 



feet eould however be erased if either the dust at the in- 
ner edge is not sufficiently heated or if the wah becomes 
permissive to molecule destroying UV radiation. In the 
following we discuss potential caveats of the model and 
how these would alter the results presented here. 

4.1. Accretion Heating 

The original model was computed assuming passive 
heating by the central star and heating due to disk accre- 
tion was treated as negligible. Accretion heating, how- 
ever, predominantly in creases the midplane temperature 
(|D'Alessio et al.l 11998). and, if significant, will increase 
the thermal desorption of molecules from grains, enhanc- 
ing the predicted effect. Details of the dynamical trans- 
port of gas and dust and their effects on the chemistry, 
however, are beyond the scope of this Letter. 

4.2. Dust Settling and Grain Coagulation 

For this model we have assumed an unsettled outer disk 
with grains of ISM abundance. In our analysis we have 
also explored a model that includes mixed grain growth 
of dust particles of up to 1 mm in size, without vertical 
settling. This lowers the UV opacity (cm^ g~^) by ~ 2 
orders of magnitude, which is a direct consequence of the 
fact that large grains contribute significantly to the mass 
but little to the UV opacity. The high model densities 
at the inner edge (pdust ^ 10~^^g cm^), however, are 
still of sufficient magnitude such that the wall remains 



optically thick to UV radiation for a model with mixed 
grain growth alone. 

Strong settling or removal of small grains can also make 
the disk more permeable to the UV radiation. We can 
ask how much dust would need to be removed/settled to 
allow the UV photons to penetrate and photodissociate 
the enhanced abundance of molecules present at the wall. 
To approximate this, the unsettled wall density is pdust ^ 
10~^^ g cm^ and the UV dust opacity is ~ 5x10^ cm^ g~^ 
at 1000 A (jWeingartner fc Drainell2001[ ). Based on the 
thermal model, the thickness of the "warm" inner edge 
is ~ 5 AU at R = 45 AU. For the UV to penetrate the 
inner 5 AU, the disk would need contain < 0.1% of the 
original amount of dust present to erase the effect seen 
here. We note however that settling will also increase the 
depth of the heating, which will in principle thicken the 
extent of the wall further. 

4.3. Dust in the Inner Gap 

An underlying feature of the disk model is the presence 
of a large inner "void" which allows the outer disk to be 
heated directly by the star. One possible explanation 
for s uch a void is grain gro wth into rocky planetesimals 
(e.g. ISkrutskie et aH I1990D . Thus one could postulate 
the existence of a small amount of undetected dust in- 
side the gap. Indeed some models infer the presence of 
some moderate mass of silicate s in the inne r disk (e.g. 
■Calvet et al...2002. . .2005; .Espaillat et al.ll2007D . This ma- 
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Fig. 4. — Simulated ALMA observations (computed with one hour integration time) of ^^CO and H2CO for a disk at 140 pc, at different 
rotational transitions with increasing J. This demonstrates the increasing contrast between the inner transition region and outer disk 
emission for higher J, along with ALMA's future capabilities at these frequencies. 



terial has the abihty to shadow the outer disk from being 
directly heated by the star, thus inhibiting the presence 
of the warm molecular interface presented here. 

First, one can ask how much dust would be required 
to cause the gap to become opaque (r ~ 10) to opti- 
cal radiation peaking at A = 0.7 /im as is the case for a 
T = 4000 K star. If we assume that the inner disk follows 
the outer disk profile (Section 12. ip . a total mass in dust 
of Mdust = 10~^ Mmooii inside 45 AU would be sufficient 
for the midplane to become optically thick to the opti- 
cal heating radiation. While this is an extremely small 
amount of dust, one must also consider the scale height 
of the material: h ~ 3.3 AU at R = 45 AU. Therefore, 
this cross-section is very small, and would not strongly 
shadow the outer disk and would not significantly al- 
ter the heating and UV irradiation except at the very 
central midplane region. Furthermor e, because the orig- 
inal model from iSauter et al.l ()2QQ9f ) included scattered 
light images, the presence of small grains at or above the 
disk scale height would significantly contribute to the 
amount of scattered light and not match observations 



for the model presented here. 

4.4. Gas in the Inner Gap 

By assuming gas and dust are co-spatial we intrinsi- 
cally assume the gap is empty of gas. Gas within the 
gap would not alter the original optical heating calcula- 
tion of the outer disk or the propagation of UV contin- 
uum photons, both of which have dust-dominated opac- 
ities. However, the presence of gas would strongly affect 
the propagation of Lyman-a as gas within the gap would 
likely be predominantly atomic. The presence of hydro- 
gen within the gap would cause a net reduction in the 
Lyman-a flux reaching the outer disk. This would ei- 
ther not change or marginally enhance the abundance of 
gas phase molecules at the truncation radius since the 
amount of photo-dissociating radiation is reduced. 

5. CONCLUSIONS 

We have investigated the chemistr y of a transiti o n disk 
employing the model structure of ISauter et al.l ()2QQ9l ) 
and have shown that a disk with a 45 AU inner void 
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depleted of dust and gas should have unique chemical 
properties. At the truncation radius, the midplane is 
both heated by optical photons and irradiated by the 
stellar UV and X-rays. The net effect is that the pre- 
viously hidden midplane is revealed as a region rich in 
gas-phase molecules, that should be readily observable 
by ALMA. For comparison, in a full untruncated disk 
model, the midplane is normally cold and shielded from 
the central star, resulting in substantial freeze-out of gas 
onto grains, preventing direct observation. While this 
analysis focuses on a specific disk model, these first re- 
sults provide insight into our ability to directly probe the 
unique chemistry of transition disks. The rich chemistry 
at the transition region will allow us to probe the physical 



conditions, such as density, temperature and kinematics. 
It is important to note that this model provides one pos- 
sibility, which can be readily tested with ALMA. We will 
be exploring a suite of disk models with a variety of ge- 
ometric frameworks both in dust and gas content in fu- 
ture work, to further explore the full parameter space of 
transition disks, with the aim of pioneering future obser- 
vational campaigns with ALMA to explore these exciting 
objects. 

We thank the anonymous referee for their helpful com- 
ments. This work was partially supported by NSF grant 
AST-1008800. J.S. acknowledges support by the DFG 
through the research group 759 "The formation of plan- 
ets: The critical first growth phase." 
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